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SECTION  1 
INTRODUCTION 


This  report  documents  the  strength  and  fracture  characteristics  of  HY-80, 
HY-100,  ana  HY-130  steels,  subjected  to  various  strains,  strain  rates,  temperatures, 
and  pressures.  The  result  of  this  characterization  is  the  generation  of  constants  for 
the  strength  and  fracture  models  in  the  EPIC-2  and  EPIC-3  codes.1,2 

The  following  sections  describe  the  test  data  and  the  analysis  of  the  data  used 
to  obtain  constants  for  both  the  strength  and  fracture  models. 
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SECTION  2 
TEST  DATA 


The  test  data  were  obtained  from  torsion  tests  at  various  strain  rates, 
t  Hopkinson  bar  tension  tests  at  various  temperatures,  and  quasi-static  tension  tests 

with  various  specimen  geometries  (notched  and  unnotched).  The  torsion  tests  were 

girformed  with  the  Biaxial  Testing  Machine  at  Southwest  Research  Institute,  the 
opkinson  bar  tests  were  performed  at  Eglin  Air  Force  Base,  and  the  quasi-static 
tension  tests  were  performed  at  Honeywell.  The  HY-80,  HY-100,  and  HY-130 
material  was  provided  by  the  Naval  Surface  Warfare  Center  (NSWC).  Physical 
properties  are  given  in  Table  1.3  Metallurgical  data  are  shown  in  Appendix  A. 

Torsion  data  are  shown  in  Figures  1  to  3.  Comparable  data  for  12  other 
materials  are  presented  elsewhere.4,5  Some  desirables  features  of  this  testing 
technique  are  that  the  state  of  stress  in  the  specimen  is  well  defined,  large  strains 
can  be  achieved  without  forming  geometric  instabilities,  and  a  wide  range  of  strain 
rates  can  be  obtained  with  the  same  testing  technique.  It  appears  that  adiabatic 
thermal  softening  is  causing  the  lower  stresses  at  the  higher  rates. 

Figures  4  to  6  show  Hopkinson  bar  test  data  at  various  temperatures.  The 
elevated  temperatures  are  obtained  by  surrounding  the  in-place  test  specimen  by  an 
oven  sucb  that  the  temperatures  are  applied  for  several  minutes  prior  to  testing. 
Although  it  is  possible  to  test  materials  to  greater  strains  than  those  shown  in 
Figures  4  to  6,  the  Hopkinson  bar  data  cannot  be  accurately  evaluated  after  necking 
begins  in  the  tensile  specimens.  Furthermore,  at  larger  strains  the  effects  of 
adiabatic  heating  can  also  complicate  the  results.  The  elevated  temperatures  show  a 
distinct  softening  effect  on  the  strength  of  the  materials. 

Photographs  of  the  fractured  Hopkinson  bar  specimens  are  shown  in 
Figures  7  to  9.  The  cross-sectional  area  of  these  surfaces  can  be  used  to  obtain 
fracture  characteristics. 

„  Quasi-static  stress-strain  data  for  both  tension  and  torsion  tests  are  shown  in 

Figures  10  to  12.  An  equivalent  tensile  flow  stress  is  obtained  from  the  torsion  data 
by  using  the  von  Mises  flow  rule,  the  tepsile  stress  is  o  =  Vjk,  and  the 
corresponding  tensile  strain  is  e  =  y/V3. 

The  stress  for  the  tension  test  data  is  based  on  the  current  area  of  the  neck,  and 
the  strain  is  defined  as  In  (AJA),  where  Ao  and  A  represent  the  initial  and  current 
areas  of  the  neck.  The  notched  tension  bars  give  a  concentration  of  hydrostatic 
tension  in  the  test  specimen.6 

This  is  clearly  shown  in  the  data  of  Figures  10  to  12.  For  a  specified  true  strain, 
the  stress  (including  hydrostatic  tension)  increases  for  the  notched  geometry.  It  can 
also  be  seen  that  the  presence  of  the  hydrostatic  tension  significantly  decreases  the 
strain  at  which  the  material  fractures. 
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TABLE  1 .  PHYSICAL  PROPERTIES  OF  HY-80,  UY-100,  AND  HY-130  STEEL 


HY-60  STEEL 

HY-100  STEEL 

HY-130  STEEL 

HARDNESS 

(ROCKWELL) 

C-21 

C-25 

C-30 

ELASTIC  PROPERTIES 

ELASTIC  MODULUS,  E 

(GPa) 

207 

207 

207 

POISSON'S  RATIO,  U 

.30 

.30 

JO 

SHEAR  MODULUS,  E/2 {1  +  v) 

{GPa) 

79 

79 

79 

BULK  MODULUS,  E/3  (1  -21) ) 

(GPa) 

172 

172 

172 

THERMAL  PROPERTIES 

DENSITY,  p 

(kO/ms) 

7746 

7746 

7685 

CONDUCTIVITY,  k 

(W/mK) 

34 

34 

27 

SPECIFIC  HEAT,  c  p 

(J/kgK) 

502 

502 

489 

DIFFUSIVITY,  k/pcp 

(m  a/s) 

.000009 

,000009 

.000007 

EXPANSION  COEF.  (VOL,),  a 

IK'1) 

.00001 1 

.000014 

.000013 

MELTING  TEMPERATURE,  T^ 

(K) 

1793 

1793 

1793 

N074S0JF 


STRESS,  T  CMP*} 
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AVERAGE  SHEAR  STRAJN,  7  -  R0  A. 


FIGURE  1.  STRESS-STRAIN  DATA  FOR  TORSION  TESTS  AT  VARIOUS  STRAIN 
RATES  FOR  II Y-80  STEEL 
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TRUE  TENSILE  STRAIN,  e=  ( 1/1 0) 


FIGURE  5.  STRESS-STRAIN  DATA  FOR  IlOPKINSON  EAR  TESTS  AT  VARIOUS 
TEMPERATURES  FOR  II Y- 100  STEEL 
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TRUE  TENSLE  STRESS,  o  (MPa) 
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TRUE  TENSILE  STRAIN,  e=  Am  (  A/A0) 


FIGURE  6.  STRESS-STRAIN  DATA  FOR  HOPKINSON  BAR  TESTS  AT 
VARIOUS  TEMPERATURES  FOR  MY-130  STEEL 
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1  mm 


FIGURE  7.  PHOTOGRAPHS  OK  Til  15  FRACTURE  PROFILES  ON  THE  HY-80 
STEEL  HOPK1NSON  BAR  SPECIMENS 
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FIGURE  9.  PHOTOGRAPHS  OF  THE  FRACTURE  PROFILES  ON  THE  IIY-130 
STEEL  HOPKINSON  HAR  SPECIMENS 
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NOTCHED  TENSION  BAR  DATA 


R  ■  2.BMM 


TEST 

SPECIMEN 


/  SMOOTH  BAR  TENSION  TEST 
t  DATA  002  S"1) 


/  EQUIVALENT  TENSILE  FLOW 
1  STRESS  FROM  SMOOTH  BAR 
T  TENSION  TEST  DATA  (BRIDGMAN) 


-  EQUIVALENT  TENSILE  FLOW 
i  STRESS  FROM  TORSION  DATA 
(VON  MISES) 


TORSION  DATA 
(^  « .01  S"1) 


HY-80  STEEL 


TRUE  PLASTIC  STRIAN 


FIGURE  10.  STRESS  STRAIN  DATA  FOR  QUASI  STATIC  TENSION  AND 
TORSION  DATA  FOR  II Y  80  STEEL 
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TRUE  STRESS  (MPa) 
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FIGURE  12.  STRESS-STRAIN  DATA  FOR  QUASI-STATIC  TENSION  AND 
TORSION  DATA  FOR  II Y -130  STEEL 
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SECTION  3 
STRENGTH  MODEL 


The  origin  of  the  strength  model,  and  the  techniques  used  to  extract  constants 
for  the  model,  are  described  in  References  7  and  8.  The  analysis  which  follows  is 
essentially  identical  to  that  described  in  these  references. 

The  model  for  the  von  Mises  flow  stress,  o,  is  expressed  as 


o  =  |A  +  Bun||l  +  Cine* ] 1 1  -  T*m| 


(1) 


where  e  is  the  equivalent  plastic  strain,  £*  =  t/i0  is  the  dimensionless  plastic  strain 
rate  for  £0=  1.0s"1  and  T*  is  the  homologous  temperature.  The  five  material  con¬ 
stants  are  A,  B,  n,  C,  m.  The  expression  in  the  first  set  of  brackets  gives  the  stress  as 
a  function  of  strain  for  t*  =  1.0  and  T*  =  0.  The  expressions  in  the  second  and  third 
sets  of  brackets  represent  the  effects  of  strain  rate  and  temperature,  respectively. 
The  basic  form  of  the  model  is  readily  adaptable  to  most  computer  codes  since  it  uses 
variables  (c,  £*,  T*)  which  are  available  in  the  codes. 

The  first  step  in  the  process  is  to  determine  the  constants  in  the  first  set  of 
brackets.  A  is  the  yield  stress  and  B  and  n  represent  the  effects  of  strain  hardening. 
Since  the  torsion  data  in  Figures  1  to  3  include  the  strain  rate  of  interest  (t*  —  1.0), 
it  is  a  straightforward  procedure  to  obtain  the  appropriate  constants  for  this  strain 
rate. 


The  same  three  constants  can  also  be  derived  from  the  quasi-tension  data  of 
Figures  10  to  12.  These  results  are  for  t*  =  0.002  and  must  therefore  be  adjusted 
for  t*  =  1.0.  This  can  readily  be  accomplished  after  the  strain  rate  constant,  C,  is 
determined. 

Figure  13  shows  the  effect  of  strain  rate  for  tension  and  torsion.  In  both  cases 
the  stress  increases  as  the  strain  rate  increases.  The  results  are  shown  for  relatively 
small  strain  values  (c  =  0.08,  y  =  0.20)  to  ensure  that  the  strain  rate  effect  is  not 
significantly  altered  by  adiabatic  thermal  softening.  Strain  rate  constants  for  both 
tension  and  torsion  can  be  determined  from  the  'least  squares’  linear  approximations 
of  the  data. 

The  effect  of  thermal  softening  is  shown  in  Figures  4  to  6.  The  thermal 
softening  fraction,  Kt  =  1  -  T*m,  is  simply  the  ratio  of  the  stress  at  elevated 
temperature  to  that  at  room  temperature.  The  thermal  softening  constant,  m,  is 
shown  to  give  an  analytic  expression  which  closely  matches  the  experimental  data. 
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At  this  point  there  are  two  sets  of  constants  (tension  and  torsion)  for  strain  and 
strain  rate  effects  (A,  B,  n,  C)  and  one  constant,  m,for  the  thermal  softening. 
Averaging  the  tension  and  torsion  constants  for  strain  rate  effects  gives  a  final  set  of 
constants  (A,  B,  n,  C)  to  be  used  in  the  strength  model.  The  resulting  stress-strain 
relationships  are  shown  in  Figures  14  to  16. 

It  can  be  seen  that  the  strength  results  for  the  three  materials  are  well 
behaved.  As  expected,  the  HY-80  steel  has  the  lowest  strength  and  the  HY-130  has 
the  highest  strength.  The  yield  stresses,  A,  increase  when  going  from  the  HY-80  to 
the  HY-130;  the  strain  rate  effects,  C,  decrease  when  going  from  the  HY-80  to  the 
HY-130;  and  the  thermal  softening  effects,  m,  are  essentially  identical  for  all  three 
materials. 


FIGURE  13.  STRESS  VERSUS  STRAIN  RATE  FOR  TENSION  AND  TORSION 
TESTS 
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TRUE  TENSILE  FLOW  STRESS,  a  (MPa) 
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TRUE  TENSILE  STRAIN,  e 


FIGURE  14.  ANALYTIC  STRESS  STRAIN  RELATIONSHIPS  FOR  ISOTHERMAL 
AND  ADIABATIC  CON  DITIONS  FOR  HY-80  STEEL 
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TRUE  TENSLE  FLOW  STRESS,  C(MPa) 
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0.0  .25  0.5  .75  1.0  1.25  1.5  1.75 

TRUE  TENSUE  STRAM.  e 


FIGURE  15.  ANALYTIC  STRESS-STRAIN  RELATIONSHIPS  FOR  ISOTHERMAL 
AND  ADIABATIC  CONDITIONS  FOR  HY-100  STEEL 


TRUE  TENSILE  FLOW  STRESS,  CfO^a) 
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TRUE  TENSILE  STRAIN,  E 

FIGURE  16.  ANALYTIC  STRESS  STRAIN  RELATIONSHIPS  FOR  ISOTHERMAL 
AND  ADIABATIC  CONDITIONS  FOR  II Y-130  STEEL 
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SECTION  4 
FRACTURE  MODEL 


The  origin  of  the  fracture  model  and  the  techniques  used  to  obtain  the  appro¬ 
priate  constants  are  presented  in  References  7  and  9.  As  was  the  case  with  the 
strength  constants,  tne  fracture  constants  are  obtained  in  a  manner  similar  to  that 
presented  in  these  references,* 

To  begin,  the  damage  to  an  element  is  defined  as 


where  Ac  is  the  incremen  t  of  equivalent  plastic  strain  which  occurs  during  an 
integration  cycle,  and  cr  is  the  equivalent  strain  to  fracture,  under  the  current 
conditions  of  strain  rate,  temperature,  pressure  and  equivalent  stress.  Fracture  is 
then  allowed  to  occur  when  D  =  1.0. 

The  general  expression  for  the  strain  at  fracture  is  given  by 


cr  =  [!),  +  I)2oxpD3o*||l  +  D41  ne*J  |1  +  D6'I"*J  (3) 

for  constant  values  of  the  variables  (o*,  i*y  T*)  and  o*  <  1.5.  The  dimensionless 
pressure-stress  ratio  is  defined  as  o*  =  om/o  where  omls  the  average  of  the  three 
normal  stresses  an; t  a  is  the  von  Mises  equivalent  stress.  The  dimensionless  strain 
rate,  i *  and  homologous  temperature,  T*,  are  identical  to  those  used  in  the  strength 
model  of  Equation  (1). 

The  five  constants  are  D1...D5.  The  expression  in  the  first  set  of  brackets 
follows  the  form  presented  by  Hancock  and  Mackenzie.10  It  essentially  says  that  the 
strain  to  fracture  decreases  as  the  hydrostatic  tension,  am  increases.  The  expression 
in  the  second  set  of  brackets  represents  the  effect  of  strain  rate,  and  that  in  the  third 
set  of  brackets  represents  the  effect  of  temperature.  For  high  values  of  hydrostatic 
tension  ( o*  >  1.5),  cfis  linearly  interpolated  between  cf  at  0*  =  1.5  and  efmin  at 
o*  =  oSpaii/o,  where  oSpa!l  and  cfllin  are  input  values  of  spall  stress  and  strain.7-9 

The  first  step  required  to  obtain  the  fracture  constants  is  to  determine  the 
effect  of  the  dimensionless  pressure-stress  ratio,  0*.  This  can  be  done  by  considering 
the  quasi-static  tension  ana  torsion  data  of  Figures  10  to  12  and  17.  For  the  torsion 
data  0*  =  0  and  for  the  unnotched  tension  data  0*  =  1/3.  For  the  notched  tension 
bar,  EPIC-2  simulations  were  performed  to  compute  0*  at  the  center  of  the  specimen. 
These  results  are  shown  in  Figure  18.  For  these  computations,  an  average  pressure 
for  each  set  of  adjacent  triangular  elements  was  used  to  eliminate  any  excessive 
stiffness  due  to  the  triangular  element  formulation.  1 1 
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FIGURE  17.  AVERAGE  FRACTURE  STRAIN  VERSUS  STRAIN  RATE  FOR 
TENSION  AND  TORSION  TESTS  FOR  HY-80, 11 Y- 1 00, 

AND  HY-130  STEEL 
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0.2 


COMPUTED  RESULTS  FOR 
NOTCHED  TENSION  BAR 


HY-80 

COMPUTED  GEOMETRY  AT  FRACTURE 

- HSKKS 


•yi twm 

■OUNDARY 


- F7JWAM 

P 

HY-80 


HY-100 


HY-130 


\ 

|  HY-130 

HY-100 


BRIDGEMAN  FACTOR  FOR 
UNNOTCHED  TENSION  BAR 


0  I _ I _ I _ I _ I - I - 1 - 

0  0.2  0.4  0.6  0.8  1.0  1.2 

EQUIVALENT  STRAIN  AT  CENTER,  E 

FIGURE  18,  PRESSURE-STRESS  RATIO  VERSUS  EQUIVALENT  STRAIN  FOR 
THE  TENSILE  SPECIMENS 
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For  the  notched  tensions  bars,  the  computed  pressure-stress  ratio  begins  at 
approximately  o *  =  0.8  and  peaks  somewhere  between  o*  =  1.3  to  1.4,  depending  on 
the  material.  After  its  peak  the  pressure-stress  ratio  gradually  decreases  until 
fracture.  It  should  be  noted  that  each  curve’s  maximum  pressure-stress  ratio  is  in  the 
same  order  as  their  respective  strengths. 

Figure  19  shows  the  equivalent  plastic  strain  to  fracture  as  a  function  of  the 
pressure-stress  ratio  for  quasi-static  conditions.  For  all  three  materials  an  approxi¬ 
mate  analytic  fit  to  the  data  was  performed  by  setting  D3  =  -0.5,  and  solving  for 
D1  and  D2. 

Figure  20  shows  the  effect  of  strain  rate  and  temperature  on  the  strain  to 
fracture.  The  fracture  strains  are  expressed  as  the  ratio  of  the  Hopkinson  bar  frac¬ 
ture  strains  divided  by  the  quasi-static  tensile  fracture  strains.  The  Hopkinson  bar 
fracture  strains  are  approximately  determined  by  measuring  the  cross-sectional  area 
of  the  post-tested  specimens  as  shown  in  Figures  7  to  9.  Unfortunately,  these  areas 
cannot  be  measured  with  a  high  degree  of  accuracy.  The  temperatures  for  each  test 
range  from  the  initial  temperature  at  the  beginning  of  the  test  to  the  computed 
temperature  at  the  strain  where  fracture  occurs. 

For  example,  in  Figure  20,  for  the  first  test  performed  on  the  HY-100  material 
the  specimen  was  initially  at  room  temperature  (T*  =  0).  If  all  the  plastic  work  is 
converted  to  adiabatic  heat,  the  computed  final  temperature  is  606K  (T*  =  .21).  The 
strain  at-which  this  specimen  fractured  was  cf  =  1.13,  at  an  average  strain  rate  of 
t  •*  550s  1.  The  quasi-static  (£  «  .002s-1),  room  temperature  (T*  =  0)  tension  test  is 
used  as  the  normalizing  denominator  for  the  vertical  axis.  For  the  HY-100  material 
this  fracture  strain  is  cf  =  1.04.  The  ratio  of  these  two  fractures  strains  is  1.09  and  is 
the  value  plotted  on  the  vertical  axis.  The  same  procedure  is  used  for  the  remaining 
three  HY-100  tests  with  the  only  difference  being  the  initial  temperature  of  the 
specimen. 

Under  normal  conditions  a  linear  least  squares  fit  is  performed  on  the  data  of 
each  material  in  Figure  20,  of  which  the  intercept  is  the  strain  rate  influence  on 
fracture  and  the  slope  is  the  temperature  influence,  The  HY-80,  HY-100,  and 
HY-130  materials  acted  differently  than  expected,  displaying  a  consistent  trend  of 
reaching  a  peak  fracture  strain  ratio  between  ,15T*  and  .35T*  and  then  becoming 
less  ductile  as  the  temperature  is  increased  (Figure  20).  Possible  explanations  for  this 
behavior  is  that  the  material  becomes  strain  rate  sensitive  at  higher  tempera¬ 
tures  as  shown  in  Figure  17,  or  it  may  be  experiencing  a  phenomena  referred  to  as 
'temper  embrittlement’  where  under  certain  tempering  processes  a  metal  becomes 
less  ductile. 

The  current  fracture  model  was  not  developed  to  handle  the  phenomena  which 
is  occurring  with  the  HY  materials.  By  applying  the  linear  least  squares  fit  to  the 
data,  one  obtains  deceiving  results,  implying  that  the  material  is  very  strain  rate 
sensitive  and  loses  ductility  as  the  temperature  increases.  A  more  realistic  approach 
to  incorporate  this  data  into  the  fracture  model  is  to  obtain  an  'average’  strain  rate 
constant  and  neglect  thermal  effects.  In  graphical  terms,  fit  the  best  horizontal  line 
to  the  data  points  of  each  material  as  shown  on  the  right  hand  side  in  Figure  20. 

Now  the  fracture  constants  can  be  determined.  The  strain  rate  constant,  D4,  is 
obtained  by  averaging  the  fracture  strain  ratios  of  Figure  20  for  each  respective 
material.  The  average  strain  rate  constant  for  all  three  is  greater  than  1.0.  This 


26 


EQUIVALENT  PLASTIC  STRAIN  AT 


FIGURE  19.  FRACTURE  STRAIN  VERSUS  PRESSURE-STRESS  RATIO  FOR 


ISOTHERMAL  QUASI-STATIC  CONDITIONS 
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ratio  indicates  that  the  strain  to  fracture  increases  slightly  as  the  strain  rate 
increases.  This  same  trend  is  verified  by  the  fracture  data  of  Figure  17,  (The 
decrease  in  fracture  strains  for  the  higher  strain'rate  torsion  data  may  be  due  to 
shear  localizations  resulting  from  adiabatic  heating.4,6)  The  constants  related  to  the 
pressure-stress  ratio  (Dl,  D2,  D3)  are  those  of  Figure  19,  adjusted  from  quasi-static 
conditions  (fi  =  .002s  i)  to  t*  =  1.0.  Finally,  the  temperature  constant,  D5,  is  set  to 
zero,  due  to  the  unconforming  data  of  Figure  20.  The  resulting  relationships  are 
shown  in  Figures  21  to  23. 

Unlike  the  strength  relationships,  which  transitioned  in  a  well-behaved 
manner  from  HY-80  to  HY-130,  the  fracture  characteristics  do  not  follow  a  well- 
behaved  pattern.  The  HY-80  steel  exhibits  the  most  ductility,  as  expected,  but  the 
HY-100  is  the  least  ductile.  Generally,  it  would  be  expected  that  the  higher  strength 
HY-130  steel  would  exhibit  less  ductility.  The  reasons  for  this  apparent  discrepancy 
are  not  clear. 

Although  no  tests  were  performed  to  obtain  spall  characteristics  of  HY-80, 
HY-100,  and  HY-130,  it  is  necessary  to  designate  a  spall  stress  and  strain  to  complete 
the  model.  Based  on  comparisons  with  materials  of  similar  yield  strengths,  the  spall 
stress  and  strain  were  estimated  and  are  shown  in  Table  2. 
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PRESSURE-STRESS  RATIO,  a*  =  a  /Tt 

IH 


FIGURE  21 .  FRACTURE  STRAIN  AS  A  FUNCTION  OF  PRESSURE-STRESS 
RATIO,  STRAIN  RATE,  AND  TEMPERATURE  FOR  UY-80 
STEEL 
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PRESSURE-STRESS  RATIO,  O*  =  o-r0 

n\ 


FIGURE  22.  FRACTURE  STRAIN  AS  A  FUNCTION  OF  PRESSURE-STRESS 
RATIO,  STRAIN  RATE,  AND  TEMPERATURE  FOR  II Y-100 
STEEL 
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HY- 130  STEEL 
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FIGURE  23.  FRACTURE  STRAIN  ASA  FUNCTION  OF  PRESSURE-STRESS 
RATIO,  STRAIN  RATE,  AN I) TEMPERATURE  FOR  II Y-I30 
STEEL 
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TABLE  2.  SUMMARY  OF  STRENGTH  AND  FRACTURE  MODEL  CONSTANTS 
FOR  HY-80,  HY-100,  AND  HY-130  STEEL 


HY-80  STEEL 

HY-100  STEEL 

HY-130  STEEL 

STRENGTH  CONSTANTS  FOR 

C-[A  +  Ben][1  +Cln  £* ] [  1  -T *"' ] 

A 

(MPa) 

872 

758 

920 

B 

(MPa) 

425 

402 

333 

n 

.36 

.26 

.39 

C 

.014 

.011 

.000 

m 

1.14 

1,13 

1.15 

FRACTURE  CONSTANTS  FOR 
e,  -[D,  4-DjEXP  ][1 +D4lnl*] 

D, 

i  +  d#t*] 

-.85 

-.61 

-.69 

Da 

2.70 

2.07 

2.39 

Dj 

-0.60 

-0.50 

•0.50 

.006 

.010 

.004 

D, 

0.0 

SPALL  STRESS,  crSPAL1 
mh 
ei 

(MPa) 

5500 

5600 

5900 

.033 

.034 

.036 

C *  -  <Jm/3 

e*-  £  /c0  FOR  E0  -1.0S’1 

T  “  H"  -  TrooW  )  f  rr  MELT  •  Tnoou ) 

N07-021 JF 
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SECTION  5 

CONSTANTS  FOR  EPIC-2/EPIC-3  DATA  LIBRARY 


A  summary  of  the  strength  and  fracture  model  constants  is  given  in  Table  2. 
Data  statements  for  the  EPIC-2/EPIC-3  material  library  (Subroutine  MATLIB)  are 
shown  in  Figure  24,  and  the  corresponding  Preprocessor  output  is  shown  in 
Figures  25  to  27. 
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c 

c  ***  HY-80  STEEL  *********** •  ••* ********************** *  TJH  -  AUG  1987 
C  **•  LIMITED  TO  U.S.  GOVERNMENT  AGENCIES  ONLY 

C  PRESSURE  -  KOHN,  AFWL-TR-69-38 ,  1969,  P.  140  (STAINLESS  STEEL  304L) 

C  STRENGTH  -  H0LMQUI3T,  NSWC  REPORT  1987 

C  FRACTURE  -  HOLMQUIST,  NSWC  REPORT  1987 

C  SPALL  -  OVERESTIMATED  FROM  DATA  IN  ArWAL-TR-81-4040 ,  1901 

DATA  MTYPE(25)/l/,IDAM(25)/0/,irAIL(25)/0/,ErAIL(25)/999./, 

1  DEN (25)/. 000725/,SPH( 25 )/4 33000./, CDUCT( 25 )/4 . 32/, 

2  ALPHA  ( 25 )/. 0000 06 3/, TEMPI ( 25)/70 ./,TROOM{ 25 )/70 ./, 

1  TMEL1(  25  )'/2768. /,  C0  (25)/.  02/,  G(  25 1/11240000./, 

4  Cl(25)/97460./,C2(25)/61640./,AN(25)/.36/,C3<25)/.014/, 

5  AK( 25)/l. 14/, C4{ 25  )/0. /,C5( 25 )/0./,C( 25 1/23776  000  ./, 

6  D( 25 )/42693000./,S( 25 1/7251 3000 ./,GRUN( 25 1/1,16/ , 

7  Q1 (25)/, 2/, Q2 (251/4./,  PMIN(25)/1 0000000. /,C9(25)/0./, 

6  D1 ( 25 )/-. 85/,D2 ( 25 )/2 . 70/, D3( 25 )/-. 50/, D4 ( 25 )/. 006/, 

9  D5 ( 25 )/0 . 00/, PFAIL( 25 1/795000 ./, EFMIN( 25 )/. 033/, C10( 25)/0 ./, 

A  (DESCM(25,J),J-1,4)  /4H*HY-,4H80  S, 4HTEEL, 4H  /, 

B  (DESCM(25,J) ,J«5,8)  /4H  ,4H  ,4HRC-2,4Hl  /,, 

C  (DESCM(  25,J),v7»9,12)/4H  ,4H  ,  4H  ,4H  /' 


C 

c  ***  HY-100  STEEL  *************************************  TJH  -  AUG  1987 
C  ***  LIMITED  TO  U.S.  GOVERNMENT  AGENCIES  ONLY 

C  PRESSURE  -  KOHN,  AFWL-TR-69-38,  1969,  P.  14U  (STAINLESS  STEEL  304L) 

C  STRENGTH  -  HOLMQUIST,  NSWC  REPORT  1907 

C  FRACTURE  -  HOLMQUIST,  NSWC  REPORT  1987 

C  SPALL  -  OVERESTIMATED  FROM  DATA  IN  ArWAL-TR-81-4040,  1981 

DATA  MTYPE ( 26 )/l/ , IDAM( 26 )/0/, I FAI  L( 26 1/0/, EFAIL( 26 )/999 ./, 

1  DENI  26 )/.  000725/, SPH( 26 )/43 3000. /,CDUCT( 26  1/4. 32/, 

2  ALPMA( 26 )/. 0000077/,  TEMPI ( 26 )/70./, TROOM(  26 )/70./, 

3  TMELT( 26 ) /2768. /,C8(  26)/. 02/, G( 26 1/11240000./, 

4  Ci(26)/109900./,C2(26)/58300./,AN(26)/.26/,C3(26)/.011/, 

5  AM ( 26 )/l . 13/,C4 ( 26  I/O ./, C5 ( 26  I/O ./, C( 26 )/2 3776000 . /, 

6  D(  26 1/42693000. /, S( 26 1/72513000. /,GRUN( 26 )/l. 16/, 

7  Q1 ( 26  )/.  2/, Q2 { 26 )/4. /, PMIN( 26  1/10000000. /,C9( 26 1/0./, 

8  Dl(26)/-.61/,D2(26)/2.07/,D3( 26 )/- . 50/,D4 ( 26 )/. 010/, 

9  D5 (  26 )/0 . 00/,PFAIL( 26 1/819  300  ./,EFMIN(26)/.034/,C10(26)/0./, 

A  ( DESCM ( 26,J),J"1,4)  /4H*HY-,4H100  , 4HSTEE , 4HL  /, 

B  (DESCM(26,J) ,J«5,0)  /4H  ,4H  ,4HRC-2,4H5  /, 

C  ( DESCM ( 26 ,J),J“9,12)/4H  ,4H  , 4H  ,4H  / 


C 

C 

C 

C 

c 

C 

c 


*  *  * 
*** 


PRESSURE  - 
STRENGTH  - 
IRACTURE  - 
SPALL 


HY-130  STEEL  *************************************  TjH  _  1987 

LIMITED  TO  U.S.  GOVERNMENT  AGENCIES  ONLY 

KOHN,  AFWL-TR-69-38 ,  1969,  P.  140  (STAINLESS  STEEL  304L) 
HOLMQUIST,  NSWC  REPORT  1987  ' 

HOLMQUIST,  NSWC  REPORT  1937 

OVERESTIMATED  FROM  DATA  IN  AFWAL-TR-81-4040 ,  1981 
DATA  MTYPE (27 )/l/ , I DAM ( 27 )/0/, I FAIL( 27 ) /0/, EFAI L ( 27 1/999 ./, 

DEN(27 )/. 000738/, SPH( 27 1/422000. /,CDUCT( 27 )/3. 46/, 

ALPHA ( 27 )/ .0000073/, TEMPI ( 27 )/70 . /, TROOM ( 27 )/70 ./, 

TMELTi 27)/2768./,C6(27)/.02/,G( 27 1/11240000 ./, 

^ii?^^3?r0,0,<':S2(27)/48300”/'AN!27)/-39/'c3(27)/*°oe/. 

^^?>/hi5/,C4(27)/0,/'c5(27  )/0./,C(  27  1/23776000./, 

D( 27)/42693000. /, S ( 27 )/72513000. /,GRUN( 27 1/1. 16/, 

Q1 ( 27 )/  .2/ ,Q2 (27)/4./, PM IN( 27 )/l 000 0000 . /,  C9 ( 27 1/0./, 

Dr  2Zy;-69/' 02(27  )/2-39/,  D3(  27  )/-.  50/,D4  (  27  1/.004/, 

D5 ( 27 )/0.00/,PFAIL(27)/8 60000. /,EFMIN(27)/. 036/, Cl0(27)/0,/, 

( DESCM ( 27  ,J)  ,  J-l ,  4  )  /4H*HY- ,  4lll 30  ,  4HSTEE  ,  4HL  /, 

( DESCM ( 27  ,  J) , J*5 , 8  )  /4H  ,4H  ,4HRC-3,4H0  / 

( DESCM ( 27 ,J) ,J«9,12)/4H  ,4H  , 4H  ,4H  / 


FIGURE  24. 


DATA  STATEMENTS  FOR  THE  EPIC-2/EIMC-3  MATERIA), 
LIBRARY  IN  ENGLISH  UNITS 
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INPUT  DATA  F OR  SOLID  MATERIAL 


MATERIAL  NUMBER  25  *HV-80  STEEL  RC-21 


MASS/THERMAL  PROPERTIES 


DENSITY 

m 

0.725000B-03 

SPECIFIC  HEAT 

m 

0 . 433000E+06 

CONDUCTIVITY 

90 

0 . 432000E+01 

VOLUME  EXPANSION  COEF 

m 

0.630000E-05 

INITIAL  TEMPERATURE 

M 

0 . 700000E+02 

ROOM  TEMPERATURE 

'm 

0.700000E+02 

MELTING  TEMPERATURE 

m 

0 .  276800E+04 

STRENGTH  PROPERTIES 

SHEAR  MODULUS 

m 

O.1124OOE+O0 

YIELD  STRESS,  Cl 

90 

0 . 9746O0E+O5 

HARDENING  COEF,  C2 

» 

0 . 616400E+05 

HARDENING  EXPONENT,  N 

M 

0.360000E+00 

STRAIN  RATE  COEF,  C3 

m 

0 . 140000E-01 

SOFTENING  EXPONENT,  M 

m 

0 . 114000E+G1 

PRESSURE  COEF,  C4 

m 

0 , OOOOOOE+OO 

MAX  STRENGTH  (OPTIONAL) 

m 

O.OOOOOOE+OO 

EQUATION  OF  STATE 

Kl 

m 

0 . 237760E+06 

K2 

m 

0 . 426930E+0B 

K3 

m 

0 . 7251 30E+08 

GRUNEISEN  COEF 

m 

0.116000E+01 

MAX  NEGATIVE  PRESSURE 

m 

0 . 100000E+08 

ARTIFICIAL  VISCOSITY 

LINEAR  COEF 

m 

0 . 200000E+00 

QUADRATIC  COEF 

Ml 

0 . 400000E+01 

HOURGLASS  COEF 

« 

0 , 200000E-01 

FRACTURE  PROPERTIES 

DAMAGE  COMPUTED 

•» 

0 

FRACTURE  ALLOWED 

« 

0 

D1 

■ 

-0.85000QE+00 

D2 

m 

0 . 270000E+01 

D3 

m 

-0 . 500000E+00 

D4 

90 

0 . 600000E-02 

D5 

« 

0, OOOOOQE+OO 

MINIMUM  FRACTURE  STRAIN 

-i 

0. 330000E-01 

SPALL  STRENGTH 

M 

0. 795000E+06 

STRAIN  FOR  TOTAL  FAILURE 

- 

0 . 999000E+03 

EXTRA  CONSTANTS 

XI  (C9)  CERAMICS  CRUSHED 

Ml 

O.OOOOOOE+OO 

X2  (CIO) 

O.OOOOOOE+OO 

PIGURK26.  KPIC-2/KP1C-3  PKKPROCKSSOR  OUTPUT  KOUTUK 

DATA  IN  KNGMSI 1  UNITS  h’OR  1! Y-8U STKKI, 
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INPUT  DATA  FOR  SOLID  MATERIAL 

MATERIAL  NUMBER  26  *HY-100  STEEL  RC-25 


MAS 3/THERMAL  PROPERTIES 


DENSITY 

M 

0.725000E-03 

SPECIFIC  HEAT 

m 

0.433000E+06 

CONDUCTIVITY 

m 

0.432000E+01 

VOLUME  EXPANSION  COEr 

m 

0 . 770000E-05 

INITIAL  TEMPERATURE 

m 

0 .  700000E+02 

ROOM  TEMPERATURE 

m 

0 .  700000E+02 

MELTING  TEMPERATURE 

m 

0.276800E+04 

STRENGTH  PROPERTIES 

SHEAR  MODULUS 

m 

0 . 112400E+08 

YIELD  STRESS,  Cl 

m 

0 . 109900E+06 

HARDENING  COEF,  C2 

m 

0 . 583000E+05 

HARDENING  EXPONENT,  N 

m 

0 .  260000E+00 

STRAIN  RATE  COEr,  C3 

m 

0.110000E-01 

SOFTENING  EXPONENT,  M 

m 

0 . 113000E+01 

PRESSURE  COEP,  C4 

m 

0 .  OOOOOOE+OO 

MAX  STRENGTH  (OPTIONAL) 

m 

o .  oqooooe+oo 

EQUATION  Or  STATE 

K1 

m 

0.237760E+08 

K2 

m 

0 . 426930E+08 

K3 

0 , 725130E+08 

GHUNEISEN  COEr 

m 

0.116000E+01 

MAX  NEGATIVE  PRESSURE 

•m 

0 . 10000QE+08 

ARTIFICIAL  VISCOSITY 

LINEAR  COEP 

m 

0 . 200000E+00 

QUADRATIC  COEr 

m 

0 . 400000E+01 

HOURGLASS  COEF 

m 

0.200000E-01 

FRACTURE  PROPERTIES 

DAMAGE  COMPUTED 

m 

0 

FRACTURE  ALLOWED 

m 

0 

D1 

m 

-0 . 61 OOOOE+OO 

D2 

■ 

0 . 207000E+01 

D3 

m 

-0.500000K+00 

D4 

m 

0 .  lOOOOOE-Ol 

D5 

tm 

0.  OOOOOOE+OO 

MINIMUM  FRACTURE  STRAIN 

m 

0 . 34Q000E-01 

SPALL  STRENGTH 

m 

0 » 819300E+06 

STRAIN  TOR  TOTAL  FAILURE 

m 

0 . 999Q00E+Q3 

EXTRA  CONSTANTS 

XI  ( C9 )  CERAMICS  CRUSHED 

m 

0. 000000E+00 

X2  (CIO) 

0. OOOOOOE+OO 

FIGURK  26.  EHC-2/BF1C-3  PREPROCESSOR  OUTPUT  FOR  TUB  MATERIAL 
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INPUT  DATA  FOR  SOLID  MATERIAL 


MATERIAL  NUMBER  27  +HY-130  STEEL  RC-30 


MASS/THERMAL  PROPERTIES 
DENSITY 
SPECIFIC  HEAT 
CONDUCTIVITY 
VOLUME  EXPANSION  COEF 
INITIAL  TEMPERATURE 
ROOM  TEMPERATURE 
MELTING  TEMPERATURE 
STRENGTH  PROPERTIES 
SHEAR  MODULUS 
YIELD  STRESS,  Cl 
HARDENING  COEF,  C2 
HARDENING  EXPONENT,  N 
STRAIN  RATE  COEF,  C3 
SOFTENING  EXPONENT,  M 
PRESSURE  COEF,  C4 
MAX  STRENGTH  (OPTIONAL) 
EQUATION  OF  STATE 
Kl 
K2 
K3 

GRUNEISEN  COEF 
MAX  NEGATIVE  PRESSURE 
ARTIFICIAL  VISCOSITY 
LINEAR  COEF 
QUADRATIC  COEF 
HOURGLASS  COEF 
FRACTURE  PROPERTIES 
DAMAGE  COMPUTED 
FRACTURE  ALLOWED 
Dl 
D2 
D3 
D4 
D5 

MINIMUM  FRACTURE  STRAIN 
SPALL  STRENGTH 
STRAIN  FOR  TOTAL  FAILURE 
EXTRA  CONSTANTS 

XI  ( C9 )  CERAMICS  CRUSHED 
X2  (CIO) 


0.738000E-03 
0 , 422000E+06 
0. 346000E+01 
0 . 730000E-05 
0. 700000  F.+ 02 
0 . 7  00O00E+02 
0 . 2768  00E+04 

0.112400E+08 
0 . 133400E+06 
0 . 483000E+05 
0 . 390000E+00 
0 . 800000E-02 
0.115000E+01 
O.OOOOOOE+OO 
O.OOOOOOE+OO 

0 . 237760E+08 
0 . 426930E+08 
0.725130E+08 
0 . I16000E+01 
0 . 1 OOOOOE+08 

0.200000E+00 
0 . 400000E+01 
0.200000E-01 


0 

0 

-0 . 690000E+00 
0.239000E+01 
-0.500000E+00 
0 . 400000E-02 
O.OOOOOOE+OO 
0. 360000E-01 
0.860000E+06 
0 . 999000E+03 

O.OOOOOOE+OO 

O.OOOOOOE+OO 


FIGURE  27.  EPIC-2/EPIC-3  PREPROCESSOR  OUTPUT  FOR  THE  MATERIAL 
DATA  IN  ENGLISH  UNITS  FOR  II Y- 130  STEEL 
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SECTION  6 
SUMMARY 


HY-80,  HY-100,  and  HY-130  steels  have  been  tested  at  various  strains,  strain 
rates,  temperatures  and  pressure.  Tb  "?  ’vaults  of  these  have  been  analyzed  to  obtain 
strength  and  fracture  constants  for  the  corresponding  models  in  the  EPIC-2/EPIC-3 
codes.  These  constants  have  been  incorporated  into  the  EPIC-2/EPIC-3  material 
data  library. 
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APPENDIX  A 

METALLURGICAL  RESULTS  OF  HY-80,  HY- 100, 
AND  HY-130  STEEL 


The  metallugrical  samples  of  each  material  were  prepared  in  cross  section  for 
examination  in  three  planes  using  m  autopolishing  technique  and  etched  with 
3  percent  nital.  The  three  orientations  examined  were:  a  plane  parallel  to  the  plate 
surface  and  two  planes  90°  to  gadh  other  viewing  the  plate  cross  section.  Various 
degrees  of  typical  inclusions  were  noted  in  the  fine  tempered  martensitic  structures. 
The  HY-130  (higher  strength  material)  exhibited  only  spot,  no  linear  inclusions. 
Tyoical  microstructures  are  illustrated:  HY-80,  Figures  A  1  to  A-5;  HY-100, 
Figures  A-6  to  A- 10;  and  HY-130,  Figures  A-l  1  to  A-14. 
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Microstructure  photomicrographs  of  IIY80  steel  In  cross  section  at  various 
magnifications,  viewing  a  plane  perpendicular  to  the  sample  length.  Views  A  &  B 
Illustrate  tempo  rod  martensite.  Etchant:  3%  nltal. 
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